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Abstract A radiation transfer model for the top-layer of a flared dust disk around 
a low mass star is presented. The disk structure, temperature distri- 
bution and spectral energy distribution are compared with the results 
of the analytical two-layer model of Chiang & Goldreich. We find that 
both, a proper iterative calculation of the main parameters and proper 
radiation transfer calculations, are important to not overestimate the 
disk structure and its emission especially at long wavelengths. 



Introduction 

Many T Tauri stars are observed with a flat spectral energy distribu- 
tion (SED) in the near infrared (e.g. Beckwith et al. 1990). This is in 
contradiction with the standard disk which is usually modeled as a flat 
blackbody. The flat spectrum in the near infrared hints to an additional 
hot dust component in the circumstellar region of young stars. Several 
models have been proposed, e.g. a dusty halo (e.g. Miroshnichenko et al. 
1999) or a dust shell beyond the outer edge of the disk (e.g. Gregorio- 
Hetem &; Hetem 2002). Both models can fit observed SEDs. However, 
it is not clear how such a configuration of the circumstellar dust should 
form. 

We are dealing here with a different location of the hot dust compo- 
nent, namely the surface of the circumstellar disk. For the surface to 
become much hotter than the interior, it must be inclined with respect 
to the infalling stellar radiation, i.e. the disk must be flared. Such a 
flaring can happen naturally when the disk is in hydrostatic equilibrium 
in z-direction with a temperature distribution that falls off more slowly 
than r _1 (Kenyon & Hartmann 1987). 
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Chiang & Goldreich (1997, hereafter CG) developped an analytical 
two-layer model. Their flared dust disk contains an optically thick 
isothermal mid-layer and a superheated surface layer. CG provide sim- 
ple analytical expressions to calculate the SEDs which are indeed flat 
in the near infrared, and their model is used by many people to fit the 
SEDs of young stellar objects. 

Our intention was to study the reliability of these analytical two-layer 
models. We developped a radiation transfer code for the upper disk part 
to calculate the parameters that determine the disk structure and tem- 
perature distribution self-consistently, and we compare our results with 
the CG model. Note that we compare our results to the original ana- 
lytical CG model. A comparison between an improved two-layer model 
with a radiation transfer model has recently been done by Dullemond & 
Natta (2003). 

1. Description of the dust disk model 

Our disk consists of an optically thick (ry ld 2> 1) mid-layer which 
is isothermal in z-direction. The mid-layer is sandwiched by two top- 
layers which are marginally optically thick at visual wavelengths but 
still optically thin at IR wavelengths, i.e. Ty° p > 1 and Tj^ p < 1. 
Consequently, the infalling stellar light is completely absorbed within the 
top-layer. The dust particles re-radiate the energy at IR wavelengths, 
so half the redistributed energy leaves the disk into space, and the other 
half penetrates the mid-layer and heats it. As a third distinct region 
we define the disk photosphere which encloses the uppermost part of the 
top-layer and whose location is defined by the parameter h (see Sect. 1.2). 

Further, we make the following assumptions: the disk is in hydrostatic 
equilibrium in z-direction, gas and dust are well mixed throughout the 
disk, and in the mid-layer gas (g) and dust (d) are in thermal equilibrium 
at the same temperature T g = = T m id- The emission and absorption 
of the gas component of the disk is ignored. 

We restrict the discussion to passive disks, which means that the only 
heating source is the star which illuminates the disk surface. In addition, 
our (gas) disk extends down to the stellar surface with no inner hole. 
The existence of such a hole would lead to a puffed-up inner rim and self- 
shadowing effects of the disk (see Dullemond et al. 2001 and Dullemond 
2002). 

1.1 Stellar illumination of the disk surface 

For a razor-thin passive disk around a star with effective temperature 
T* and radius R* the monochromatic flux entering the disk perpendicular 
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Figure 1 Sketch of a flared 
disk. The flaring term of 
the grazing angle is indi- 
cated. The inner radius of 
the dust disk is ro- 



to the surface at distance r from the star is 
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where the star is assumed to emit a Planck spectrum. This flux can be 
parametrized in terms of the solid angle 0, under which the star is seen 
from a disk surface element at distance r and the so-called grazing angle 
a gr , i.e. the mean angle of incidence of the stellar flux 



The grazing angle of a razor-thin disk is therefore 
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which for large distances reduces to the handy formula used by CG 
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This relation also holds for a wedge-shaped disk as long as its opening 
angle is small. 

In a flared disk the surface is curved and the grazing angle increases 
with distance from the star (see the the solid lines in Fig. 1). There- 
fore, the flared disk intercepts more stellar light and is heated more 
efficiently. The grazing angle of the flared disk is by the amount a^ rc 
greater than for a flat disk (Fig. 1), and the total grazing angle becomes 
simply a gr = a T g f or + re 
can be computed from 

dh h d / 
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With this new grazing angle, which is a function of the location h of 
the photosphere, we can calculate the flux penetrating the flared disk in 
z-direction. 
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1.2 The location of the disk photosphere 

The photosphere encloses the uppermost part of the top-layer. Its 
onset is described by the parameter h, which is defined as the height z 
above the mid-plane where the visual optical depth along the direction 
of the inf ailing stellar light, i.e. along the grazing angle, equals 1. This 
leads to a vertical optical depth of 

Ty = Sin Q gr . (6) 

Since the disk is assumed to be in hydrostatic equilibrium, Ty can be 
calculated 

oo 
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where po is the density in the mid-plane and H is the scale height of the 
disk given by 

H = J-¥±- r 3 / 2 (8) 

Here, M* and p are the stellar mass and the mean molecular weight. 

From equaling Eqs. (6) and (7) h can be determined, but it is a func- 
tion of temperature and grazing angle. 

1.3 The temperature of the mid- layer 

The downwards and upwards directed fluxes, and F\ are in equi- 
librium throughout the passive disk. We can calculate these two fluxes 
explicitely at the boundary between the isothermal mid-layer and the 
top-layer. We assume that half the incident stellar flux is re-radiated 
into space and the other half penetrates the mid-layer, leading to 

Fi = \f ± = \j " <x ff SlB v {T*)dv (9) 

The flux leaving the mid-layer in upward direction is given by 

F^ = 2ir J J B u (T mid ) (l - e- r ^) p dp dv (10) 

where T m id is the isothermal temperature of the mid-layer, and we set 
p = cos 9 with 9 as the angle measured from the z-axis. The visual 
optical depth of the mid-layer in vertical direction can be written in the 
form 
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with the visual optical depth at the inner edge, Ty ld (ro), and the ex- 
ponent s as free parameters. The temperature of the mid-layer, T mi( i, 
follows from equaling Eqs. (10) and (9). This temperature is, however, 
a function of the grazing angle. 

We have now three important parameters, a gr , h, and T mi ^, but none 
of them can be computed independently, instead they must be calculated 
iteratively. 

1.4 Radiation transfer within the top- layer 

We have to specify the visual optical depth in the top-layer in z- 
direction, Ty° p . It should be high enough so that first, the infalling 
stellar radiation is completely absorbed and second, the dust grains at 
the bottom of the top-layer reach a temperature close to that of the 
mid-layer in order to guarantee a smooth transition. On the other hand, 
Ty° p must be small enough to allow the dust emission, which occurs at 
infrared wavelengths, to escape the top-layer. 

The radiation transfer equation of a plane-parallel slab, 

/„(//, r„) = / e T ^ - / SrWe-^-V" - (12) 

with the source function S u , the incident intensity Iq and fj, = cos 6 is 
split into up-streams, /+, that penetrate the top-layer from the mid- 
layer, and down-streams, I~, that cross the top-layer starting from 
the surface. The incident intensity is either the stellar radiation for 
downwards directed streams (but only for angles 6 under which the 
star can be seen, else it is zero), or the emission of the mid-layer, 

£>j,(T m id)(l — e -7 "™' for the upwards directed streams. 
With the help of the Feautrier parameters 

U„ = \{lt + K) and Vu = l(I+-I-) (13) 
the mean intensity J„ becomes 

Jv = Ju{j u ) = J U U (H, T u ) dfJL (14) 

which is needed to calculate the source function and the emission of the 
grains. The source function itself determines the up- and down-streams 
of the intensity; the radiation field calculation must therefore be iterated. 

Finally, the temperature of the grains at each location in the top-layer 
follows from balancing their absorption from the surrounding radiation 
field, Jjy, and their emission at their equilibrium temperature, 



6 




10 1 1 10 1 10 2 

r [AU] 



ZD 

^ 100 











/ - 












CG / ' - 






/ 






/ - 






/ 






/ 






^-^ow model \ 









100 200 

r [AU] 



Figure 2. Drop of the parameter 
h/H (top, see text), and the height h 
of the photosphere (bottom). 
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Figure 3. Temperature of the mid- 
layer (top), and comparison of the 
top- and mid-layer temperature at the 
boundary (bottom). 



2. Results 

Since we want to compare our results with those of CG we took the 
same parameters as they did: for the star, T* = 4000 K, L* = 1.44 L , 
R t = 2.5 R , M* = 0.5 M ; for the disk, r in ~ 0.07 AU, r out = 270 AU, 
Ty ld (r) ~ 4 • 10 5 r^u' 5 ; for the dust grains, spheres of radius a = 0.1 fim 
with mass density pd = 2gcm~ 3 and absorption efficiency Q = 1 for 
A < 2ira and Q ~ 2na/X for A > 2na. The latter implies an absorption 
coefficient Ky ~ 4 x 10 4 cm 2 per gram of dust. Scattering is absent and 
we use Ty° p = 3. The disk is seen face-on. 

We first compare our height h of the photosphere (Fig. 2) with the 
CG results. The ratio h/r which determines a gr can be written in terms 
of the scale height H 

h ~±Z (15) 
r H r 

According to CG the ratio h/H drops from 5 at 3AU to 4 at 100 AU 
but they use h/H = 4 throughout their calculations which leads to an 
extremely thick flared disk with h(r = 270 AU) = 270 AU. Our calcu- 
lations show, however, that h/H drops from about 6 at ro to < 2.8 at 
the outer edge, and our disk stays much flatter at large distances. In 
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Figure 4- Temperature distribution 
within the top-layer at r (top) and at 
r = 250 AU (bottom). The star is the 
CG value, square and triangle are our 
and the CG mid-layer value. 



Figure 5. Comparison between CG 
and our results for the top- layer (top), 
mid-layer (middle) and total SED 
(bottom). 



a recent work Chiang et al. (2001) corrected their old results for the 
disk thickness and their corrected photosphere almost agrees with our 
results. 

The temperature in the mid-layer, (Fig. 3), is not so different 

in the two models. Our mid-layer is slightly hotter especially in the 
inner parts (r < 20 AU) but comparable with CG for larger distances. 
Our temperature calculations in the top-layer result in a temperature 
gradient between the very hot surface and the bottom of the top-layer. 
In Fig. 4 we compare our temperature distribution with the CG surface 
temperature as well as with our and the CG mid-layer temperatures 
at two different distances. We cannot reproduce the very hot surface 
temperature of CG, but on the bottom of the top-layer (at Ty° p = 3) 
the temperature approaches everywhere the value of the mid-layer (see 
bottom panel of Fig. 3). 

The differences in the disk and temperature structure between the 
CG model and our more detailed radiation transfer calculations lead of 
course also to differences in the SED (Fig. 5). The emission from the CG 
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surface layer which is much hotter than our top-layer is shifted to shorter 
wavelengths. The emission of the CG mid-layer comes from a narrower 
wavelength region than in our model and is too low. The total SED of 
CG therefore deviates from our results especially in the long wavelength 
region (A > 30/xm). 

3. Conclusions 

A radiation ransfer model for the upper disk part of a flared disk is 
presented and the results are compared with the analytical two-layer 
model of CG. We show that the major parameters that determine the 
structure of the disk, a gr , h, T m ^, must be calculated iteratively. In addi- 
tion we perform detailed radiation transfer calculations in the top-layer 
to find the temperature structure in the upper disk part. At the bottom 
of the top-layer the temperature agrees very good with the values of the 
mid-layer. The very hot surface temperature found by CG could not be 
reproduced, and our total SED deviates from the one of CG especially in 
the long wavelength regime. Although the two-layer model gives handy 
formulae to calculate the SED of a flared dust disk, the results should 
only be used as a 'first guess' for the structure and temperature distribu- 
tion of the disk and radiation transfer calculations should be performed 
for a better characterization of the circumstellar dust disk. 
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